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The atomic electric dipole moment (EDM) currently searched by a number of experimental groups 
requires that both parity and time-reversal invariance be violated. According to current theoretical 
understanding, the EDM is induced by the nuclear Schiff moment. The enhancement of the Schiff 
moment by the combination of static quadrupole and octupole deformation was predicted earlier. 
Here we study a further idea of the possible enhancement in the absence of static deformation but 
qq \ in a nuclear system with soft collective vibrations of two types. Both analytical approximation and 

, numerical solution of the simplified problem confirm the presence of the enhancement. We discuss 

related aspects of nuclear structure which should be studied beyond mean-field and random phase 
' approximations. 
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Efforts of many experimental groups are directed to the discovery of the electric dipole moment (EDM) of atoms. 
The existence of the EDM in a stationary state of a finite quantum system would manifest the simultaneous violation 
of parity {V-) and time-reversal (T-) invariance [lj. The best current experimental limits @, Q are found for 129 Xe 
and 199 Hg. Further progress in this field can be connected with the better understanding of nuclear many-body 
mechanisms which could substantially enhance the effect. This problem is also of significant interest for theory of 
^ ■ collective motion in nuclei and other mesoscopic systems. 

The atomic dipole moment is induced by the weak interactions of atomic electrons with the nucleus. The dipole 
moment of the nucleus is almost completely screened by the redistribution of atomic electrons in the applied external 
electric field (the famous Schiff theorem Instead of the nuclear dipole moment, one needs to consider the vector 
operator of the next order, the so-called Schiff moment 
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The correctness of this conclusion was confirmed in the recent theoretical discussion that lead to the reexamination 
of the Schiff theorem 0, 0] . The similar operator is known as generating the isoscalar giant dipole resonance in nuclei 
since the isoscalar dipole moment creates only a spurious mode of the center-of-mass displacement. 

In the presence of V, T- violating weak interaction, the non-zero expectation value of the Schiff moment, as of any 
polar time-even vector, becomes possible in the ground state of a nucleus with non-zero angular momentum I, 



(s) = ((s ■ i)) jjY+Ij' < 2 > 

In the next step, this vector induces the EDM of the atom that is enhanced in heavy atoms. Typically, the Schiff 
moment has a single- particle order of magnitude value as confirmed by various calculations of nuclear structure, see 
for example [g, [9|, LLQi |ll| ; the core polarization effects can increase the result by a factor of order 2. 

It is however known that there exist promising collective mechanisms which must be explored in the search of 
the many-body enhancement of the Schiff moment. Parity non-conservation in scattering of polarized neutrons off 
heavy nuclei is enhanced by orders of magnitude [12j | because of extreme proximity of p- and s-wave resonances and 
the chaotic nature of their extremely complicated wave functions. In contrast, here we are interested in coherent 
mechanisms which could show up in the properties of the ground states and low-lying excitations. Essentially we 
need the collective states of opposite parity in the vicinity of the ground state which can be effectively admixed to 
the ground state by the weak interaction. The main hope here is associated with coexistence of quadrupole and 
octupole collective degrees of freedom. Indeed, it was shown that the intrinsic Schiff moment (in the body-fixed 
frame), that exists without violation of fundamental symmetries, can be enhanced by 2-3 orders of magnitude in 
nuclei with simultaneous quadrupole and octupole static deformation fl3l . Ihll . [l5| . The mixing by the weak interaction 
of states with certain angular momentum in the laboratory (space-fixed) frame will be particularly enhanced for parity 
doublets [HI [I?], HaV Those are the states with same spin and opposite parity which differ by the "right" and "left" 
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orientation of the asymmetric configuration that is characteristic for the deformation without reflection symmetry in 
the equatorial plane. 

The requirements of static octupole deformation can be relaxed since almost the same effect can be reached [lH [l9| 
when static quadrupole deformation is combined with a soft octupole vibrational mode. Similarly to the case of static 
octupole deformation, where the enhancement effect is proportional to the square /3f of the octupole deformation 
parameter, here the analogous result contains the dynamic mean square value (/3|) that has the same order of 
magnitude being proportional to the inverse soft mode frequency, I/W3. This significantly widens the set of nuclear 
candidates where one can expect the collective enhancement of the Schiff moment. 

There is a theoretical possibility that the enhancement may occur even in the absence of static quadrupole defor- 
mation due to the combination of quadrupole and octupole soft modes [20I [2lj . This would broaden even further 
the possibilities for the experimental search of the EDM adding, for example, such nuclei as light isotopes of radium 
and radon, where the presence of both soft modes with strong dipole transitions connecting corresponding vibrational 
bands is well known [22|]. Although phenomenological estimates [2(| supported this idea, the detailed microscopic 
calculation [2l[ based on the random phase approximation (RPA) and quasiparticle-phonon coupling gave essentially 
negative results. The effect of strong enhancement was found only in the unphysical limit of very low vibrational 
frequencies loi and W3. However, the applicability of the RPA to the situation of strong interaction of low- frequency 
collective modes is questionable and that result might be an artifact of the inappropriate approximation. 

In the current work we present a simplified model, where an unpaired particle interacts with two coupled soft 
collective modes. This interaction is effectively strong and creates a coherent state (condensate) of phonons although 
the even core is still spherical. Such a situation cannot be adequately treated in the RPA framework. We consider 
both the approximate analytical approach and exact numerical solution. The results clearly demonstrate the existence 
of a parameter region, where we see the strong enhancement of the nuclear Schiff moment. On the other hand, the 
whole problem is quite interesting irrespectively of the symmetry violations, for the study of nuclear structure in the 
situation of shape instability due to the presence of soft modes. A more detailed review of the entire problem of the 
Schiff moment and the previous search for its collective enhancement can be found in 23] . 



II. MICROSCOPIC JUSTIFICATION OF ENHANCEMENT 
A. Spontaneous symmetry breaking in an even nucleus with soft modes 



We consider a generic model of a "soft" nucleus that has well developed vibrational modes still keeping, on average, 
a spherical shape. Although the non-zero Schiff moment requires non-zero nuclear spin and therefore an odd-A 
nucleus, we start with the neighboring even-even nuclei. Here we assume that the collective modes with spin-parity 
characteristics 2 + and 3~ have low excitation energies L02 and W3, well below the gap energy 2 A that separates the 
ground state + from the threshold of two-quasiparticle states. This situation, for example, can be found in the 
chain of even-even xenon isotopes and, most clearly, in light isotopes of radon and radium [22], where we see 
soft quadrupole and octupole collective modes evidently correlated in their energetics. As the energy spectrum does 
not display a rotational pattern, we can assume that these nuclei are still spherical. However, the spectra do not 
show usual phonon multiplets either. In radium and radon isotopes, strong dipole transitions are present between 
the corresponding members of long quasivibrational bands of positive and negative parity. This situation seems to be 
favorable for dipole correlations, and consequently for the Schiff moment. 

The correlation between quadrupole and octupole vibrations is an anharmonic effect that can be theoretically 
described only if one goes beyond the random phase approximation (RPA). In the following part of the article, Sec. 
Ill, we show a numerically solved two-level model that starts from single-particle mean field levels and an inter-particle 
interaction; then both phonons and anharmonic effects result from the exact diagonalization. In the present Section 
we use a simpler analytical approximation of bosonic phonons interacting with unpaired particles. The model is a 
generalization of the approach suggested in Ref. [25[ for a global analysis of quadrupole-octupole correlations. This 
idea was later used for the explanation of such correlations seen in the experimental data for xenon isotopes [2~i| and 
in the schematic arguments applied to the ehnancement of the Schiff moment in Ref. (2lj . The correlations between 
quadrupole and octupole modes can lead to new observable effects and phase transitions also in nuclei with A > 220 
as it was discussed recently [HI, H?} ■ 

The model includes two groups of spherical single-particle levels with spins {j} and {j 1 } of positive and negative 
parity, respectively. The elementary bosonic quanta, phonons, are coherent superpositions of many two-quasiparticle 
excitations. The quadrupole phonons contain (jij^)^^ and (jiJ^p combinations allowed by the angular momentum 
coupling, while the octupole phonons are built of the pairs (Jij^a^, based on the orbitals of opposite parity. The 
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effective Hamiltonian of low-lying states in the even nucleus is given by 

H = H 2 +H 3 + H 23 , (3) 

where H 2 and H 3 describe quadrupole and octupole collective modes (in principle including their anharmonicity) . 
The interaction term, 



E[(/ f /)2,4+i>4 ( 4 ) 



#23 = X 

is responsible for the mode- mode interaction with a; as a coupling constant. The construction of ^ is the simplest 
one that accounts for parity and angular momentum conservation. 

The non-linear interaction (01 leads to the effective deformation of the excited states when both modes are simulta- 
neously present. Indeed, the low-lying quadrupole and octupole modes support each other and generate the coherent 
state with the non-vanishing expectation value of the quadrupole moment according to 

(^>=--((/ t /) 2 , I ), (5) 

Ul 2 

as follows from the operator equation of motion 

[d^H] =uJ2d^ + x(pf)^. (6) 

for the quadrupole phonon operator d^. For simplicity, here we neglect the possible quadrupole anharmonicity 
that would only substitute the unperturbed quadrupole frequency w 2 by the effective curvature of the ground state 
quadrupole potential [28[. 

The result ((5|) means spontaneous violation of rotational symmetry by the soft quadrupole mode that is determined 
by the octupole phonon which selects the intrinsic axis (we assume axial symmetry). The direction of the effective 
self-consistent deformation is arbitrary, and, to restore the symmetry and appropriate quantum numbers of total 
nuclear spin in the space-fixed coordinate frame, we accept that the orientation of the deformation is given by a 
spherical function of corresponding rank, 

<d M )=£ 2 ^Y 2 ;(n), (U)=f^Y&(n). (7) 

Here n is the unit vector of the symmetry axis considered as a variable in the collective space. A similar operator 
approach was used long ago in the derivation of the nuclear moment of inertia without applying a cranking model 

The number N$ = = f'f of octupole phonons is conserved, -/V3 = 1 in the lowest 3~ state, and f* is the 

operator generating the octupole vibrational mode in the body-fixed frame defined by the orientation n. Then eq. 
([5]) equates the n-dependence and, with the ansatz ([7]), provides the effective quadrupole deformation parameter 62, 
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The 1/w-dependence of the effective deformation parameters is the characteristic feature of the whole consideration 
based on the assumption of soft collective modes. 

The equation of motion for the octupole mode, given by the commutator and the effective quadrupole 

parameter (jHJ), is linear. It relates the excitation energy £3 of the octupole phonon with the corresponding unperturbed 
energy u> 3 and the quadrupole condensate ©. Collecting again the terms expressing the angular dependence, we obtain 

8 x 2 

E 3 =^--—. 9 

21 UJ 2 

This simple regularity first discussed in Ref. [25| in the global review of octupole vibrations provides a clear correlation 
between the two modes. Recent measurements for the long chain of even-even xenon isotopes [H| show precisely such 
a correlation, with a rather large magnitude for the parameter x that exceeds the expectations for the anharmonic 
mode-mode coupling based on the standard RPA estimates. We assume that the Schiff vector S has a non-zero 
reduced matrix element S° between the one-phonon states of the two modes, |2 + ) and |3 _ ). 
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B. Phonon condensate in an odd nucleus 

Now we consider the situation in the adjacent odd- A nucleus, where the presence of the unpaired particle in the 
ground state breaks symmetry for both soft modes. The state of the odd particle is defined by the density matrix 



Pjim 1 ;j 2 m 2 ~ a j 2 m 2 a 3l"n 



(10) 



in terms of the operators of creation, aj m , and annihilation, a,j m , of the particle. The expectation value (p^miy^mj) 
in the ground state with spontaneously broken symmetry can be written as 



LA 



(11) 



Here the even-L parts come from the pairs of levels (ji, J2) or (j'1,3'2) of the same parity, whereas the odd-L ones 
correspond to the combinations (ji, j' 2 ) and (31,32) of single-particle levels of opposite parity. 

In the odd nucleus, the unpaired particle interacts with the soft modes, and the phonon Hamiltonian should be 
supplemented with 



del 



h + £T<+> + ff'W + 



(12) 



where h contains unperturbed spherical single-particle energies 6j and , or quasiparticle energies if pairing is 
included. The positive-parity parts are given by 



H (+) -- V x {+) o 



jij 2 m 1 m 2 [i 



,/'-/j-'»a [ Jl 2 32 

mi — ji —m,2 



(13) 



and H ,ty+ \ where all positive parity levels 3 are substituted by the negative parity levels j' . The Hermiticity conditions 
for the Hamiltonian (jlip and its negative-parity analog read 



The interaction with the octupole phonon that changes parity of the particle is, analogously, 

^ = — yjlj 2 Pjim 1 ;j' 2 m 2 
jlj' 2 m 1 m 2 H 



(14) 
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mi — /i 
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h.c. 



(15) 



Here we take the coordinate part of the phonon field, for example, (dj, + (— )^d_ (tl J , since it is proportional to 1/ 

in contrast to the momentum part (with sign minus between the phonon creation and annihilation operators) that is 
small, cx y/uj, for a soft mode. 

With the directions specified by the unpaired particle, the soft modes acquire the condensate components in the 
ground state of the odd nucleus. In the same way as shown by eq. ([7]) we find 



(d„> = ftl^n) = — 



3132 
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(16) 
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(17) 
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C. Density matrix of an unpaired particle 



The presence of phonon condensates self-consistently creates effective deformation of the mean field for an unpaired 
particle. This is revealed in the appearance of deformed components pl^q of the density matrix pip . The equations 
of motion for the single-particle operators can be written as 



[a jm ,H] = e jajm - £ 4 H J(4 + (-)^)(-)" +i - m (^ 1 4 i) 



- E yf4ti + (-rf-»)(-r +j - m (i 4 i)a jim , (is) 

A similar equation is valid for the operators aj/ m related to the negative parity levels with obvious substitutions 
x (+) -m/M and //, , • //,,;. 

As now the odd particle is moving in the mean field of effective deformation, its angular momentum j is not conserved 
anymore but, owing to the assumed axial symmetry, the angular momentum projection k onto the symmetry axis still 
characterizes the single-particle levels. The appropriate transformation for the single-particle operators is therefore 

a jm = ^ 1 DmK C jK> (1^) 

K 

where we use the matrix elements D J mK of finite rotations. Such angular functions, after straightforward algebra, 
appear in all terms of equations of motion (fTSf justifying the ansatz ([19]) . The new operators Cj K satisfy, instead of 
(HHJ) , the set of equations 

it ^ ' 



- 2 ^v£E^H ji - K (ioO c ^ (2o) 

and analogously for the operators cy K for the levels originally of negative parity. In the intrinsic frame parity is mixed 
by the simultaneous presence of quadrupole and octupole deformation [we used in eq. (f2TJ)) the static values p^|) and 
P7p]. These coupled equations provide deformed single-particle states and their energy spectrum. For the states with 
one unpaired particle, the transformation (|19|) is normalized as 

E(i^) = E(i%) = i- (2i) 

jm jK 

When the linear set of equations ([20]) is solved and the Nilsson-type single-particle orbitals with no definite parity 
are found, we can return to the density matrix (jlip and make the consideration self-consistent. From eq. (jlip we 
obtain 

p L {hj 2 ) = y/^(2L+l)(-) L E(-) J2 " K ( L 2 K 3 l o ) ( C L C ^- ( 22 ) 
A similar expression is valid for Phij'xj'-i) and for the mixed-parity components of the density matrix PlOiJ^)- 



D. Schiff moment 



The operator S\ v of the Schiff moment in the even nucleus has a reduced matrix element S° = (2||Si||3) between 
the low-lying 2 + and 3~ states. As mentioned earlier, the dipole transitions between the states of the quadrupole 
and octupole bands are empirically known to be enhanced in nuclei of our interest, such as light radium and radon 
isotopes [22j . The collective contribution to this operator can be written in terms of our phonon variables as 

s lv = s° ER 1 ^ ( l v 4 I ) (<Jy + (-^"'/V-*) ■ (23) 
tip.' 
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With the ground state expectation values of the effective deformation parameters in the odd nucleus, eqs. (fT6|) and 
(jTTJ) , this gives a rotational operator 

(24) 

enhanced by small collective frequencies. As a result, we reduce the whole problem to that of the "particle 4- rotor" 
type [LH . |30| , where the static deformation is substituted by the effective deformation coming from the soft quadrupole 
and octupole modes of the spherical even core. 

The observable Schiff moment in the laboratory frame can come only from the explicitly acting V- and 7"- violating 
weak interaction W that creates an admixture of the states \n) of opposite parity to the ground state |0), 

{ ' ^ E -E n [ ' 

n 

The states \n) must have the same spin / ^ as the ground state of the odd nucleus since W is a rotational scalar. 

E. Simple model 

The structure of equations becomes clear in the simple model case of two levels j and j' of opposite parity but the 
same magnitude j — j' . The single-particle spherical energy levels are e and e', while the corresponding single-particle 
operators in the intrinsic frame are c K and c' K , where k is the projection of the single-particle angular momentum onto 
the symmetry axis. The analysis essentially follows the consideration [14j developed for a statically deformed nucleus 
with no symmetry with respect to the reflection in the equatorial plane. 

The effectively deformed core is presented by rotational states (D^ [IQ )* = \J 'Air/2L + 1 Ylm' , where parity is simply 

(— ) L . The particle orbitals are defined by the operators cj. for positive parity and c'J for negative parity multiplied 
by the corresponding rotational functions D J mK which are coupled with the core Z) L -functions to the total angular 
momentum I and laboratory projection M of the nucleus (the body- fixed projection is K = k). Due to the presence of 
octupole deformation, the intrinsic particle states | k) do not have certain parity. Their wave functions can be written 
as 

l«) = (64 + ^)10), tl + r& = l, (26) 
and the amplitudes £ K and rj K satisfy the secular equations with split energies E K , 

[E K - e + b 2 (K)}^ + b 3 { K ) VK =0, (27) 

h(K)£ K +[E K -e' + b 2 (K)]r] K = 0. (28) 
Here the effective deformation amplitudes 62(f) an d 63 (n) are given by 

62(f) = \[^P2 X 2 , 63(f) = yj/?3 X 3 . (29) 

The energy spectrum of deformed orbitals is obviously 

El ±] = e -±f - 62(f) ± \ ^(6 - £ ') 2 + 46|(f). (30) 
The deformation parameters are found from eqs. (|lfc>[) and (jTTJ) , 

&L = / or — PL: (31) 
(2L + 1)u)l 

where xi — + a;'-,,-/ ^ and X3 = 2yjji . The density matrix (|22p is now determined by the quantum numbers k of 
the actually occupied pair of orbitals, 

Pl {k) = (-)VM2I + l)(-r ( 1 K 1 ) z l(k), (32) 

where 

22(f) = £ + r)l = 1, z 3 (k) = 2C K ry K , (33) 

and in the degenerate limit (e = e') we have |£ K | = \r) K \ = l/y2, \z 3 (k)\ = 1. Due to time-reversal invariance, the 
same energy splitting (|30p occurs for the orbitals with the opposite sign of projection k. 
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F. Parity doublets 

The single-particle intrinsic states \k) and | — k) found above, eq. (I26p . should be combined into intrinsic orbitals 
of certain parity, 

\rl>* ; ±) = -j=[\K)±(-) i -' t \-*j\- (34) 

Since in the model each state \k) or | — k) is a combination of orbitals based on mixed levels j and j', in fact we 
have a quartet of relevant single-particle states or two parity doublets with the same |k|. If the matrix element of the 
weak interaction between the single-particle spherical states \ jm) and \ j'm) equals w, its analog in the intrinsic frame 
connecting the partners of the doublets (|34p is wt; K r) K . As it is clear from the symmetry arguments and seen from the 
set of equations (|27l28p , this product is proportional to the octupole effective deformation (3% and the component p3 
of the single-particle density matrix. However, this quantity becomes a constant in the exceptional case of degenerate 
orbitals, eq. (|32[) . Then the maximum value of the single-particle contribution of the weak interaction is, in this 
model, just \w\/2. 

The Schiff moment operator (|2"4")) is responsible for the matrix element between two rotational states with the same 
spin I but opposite parities. This geometric matrix element gives, as follows also from general arguments [Bl [lil [20|. 

(IMk, U\S 10 \IMk, -II) = ~ S°p 2 p 3 ■ (35) 

V7T 1(1+1) 

The final answer (|25|) can be obtained multiplying the result l[35|) by w^ K ri K / AEi K , where the energy denominator 
is given by the splitting of rotational partners of opposite parity. The collective enhancement for small frequencies 
of the even core is explicitly present in eq. p5[) . The energy splitting AEi K is determined mainly by the distance 
between the headbands of opposite parity that is of the order of the difference of phonon frequencies Auj — \ui2 — ^3]. 
The effective moments of inertia of bands in the odd nucleus built on the parity doublets are expected to be close in 
heavy nuclei. For soft vibrations, Aw gives another enhancement factor. 

III. NUMERICAL DIAGONALIZATION OF THE MODEL 

Here we show the results of the exact numerical solution of the two-level model based not on the assumption 
of preexisting soft phonons but on the full account of inter-particle interaction, including the primary interaction 
responsible for the existence of soft modes. We will be looking for the region of the parameter space that is favorable 
for the enhancement effects. 

A. Nuclear Hamiltonian 

We consider, as above in Sec. 2.5, the single-particle space of two spherical levels with the same j = j' and opposite 
parity; their spherical energies are taken as €j = and eji — e. The dynamics in the particle-particle channel can be 
expressed in terms of the pair operators from the same j-level, 

= —7= ^jmjm' a jm'a.j m , PlM = ^jmjm' a j'm'Clj'm, (36) 

mm' mm' 

where only even total momenta, L = 0, 2j — 1, are allowed, and the pair operator changing parity of the state, 

RLM = 2_j < ~'fm'jm a jma.j>m' 1 (37) 
mm' 

where all values of the pair spin, L = 0, 1, ...,2j, are possible. The most general two-body Hamiltonian, 

H lBt = H jj +H j , j ,+H jj ,+V, (38) 
includes the interactions inside each j-level, 

Hn = E h L P f LM PLM, H fj , - Yl ^lP'lmP'lm, (39) 

LM LM 



8 



the pair transfer between the levels, 

Hjf = ]T 9 L (p ] lm P'lm + P'lm^Plm) , 

LM 

and, finally, scattering of the mixed pairs, 

V = lE VlRImRlm. 

LM 



(40) 



(41) 



We specify the model choosing only the interactions associated with the collective dynamics: pairing part G > 0, 

h = K = g Q = -G- (42) 

quadrupole-quadrupole interaction in the particle- hole channel that, after the angular momentum recoupling, corre- 
sponds to 



fiL = h' L = 5x 2 

and octupole-octupole interaction that similarly gives rise to 

V L = 7x 3 



3 3 2 \ . 
3 3 L 



i i 3 

3 3 Lf 



(43) 



(44) 



Thus, the model Hamiltonian is characterized by four parameters, e, G, x 2 and x 3 . Attractive forces correspond to 
positive G, x 2 and x 3 . 

Also useful are the operators in the particle-hole channel, 



h A j 2 \ t 



-mi fj. m 2 I a nm 1 a 32m 2 , 



for different angular momenta A and all combinations of j\ and j 2 ■ These operators have a symmetry property 
In these terms, our multipolc-multipole interactions are 



H (2) = _ 



hx 2 



(45) 



(46) 



(47) 



and 



ff( 3) = _Z|3 [Ql(fj)Q 3ll (j'j) + Ql(jj')Q 3 ,(jf) 



The Schiff moment operator is given by 
while the weak interaction, 



W = ^/2j + lw Qoo(i'j) + Qoo(jj') 



(48) 



(49) 



(50) 



is characterized by one scalar w for the transition between the states \jm) and \j'm). As we are interested in the 
trends of the expectation values of the Schiff moment in a function of nuclear forces, the values of the coefficients s 
and w are irrelevant. Below we set w = 1 and (j'm = j\Sfj, = o\jm = j) = \fj] {j + 1). 
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Figure 1: The spectra of the 6-particle system for e = 0.2 and G = 0.2. For the left four spectra, a fixed value of the octupole 
strength, xz = 1.0, and a variable quadrupole interaction parameter, X2—0.5, 1, 1.5 and 2, are used. The most right spectrum 



corresponds to the pairing Hamiltonian, xi = xs = 0, where seniorities (sj 
are shown in red. 



Sji) are also marked. The negative parity states 



B. Numerical results 



Fig. [T] shows the results for the energy spectrum obtained in the exact diagonalization of the strong Hamiltonian 
for 6 particles on two levels with j = j' = 17/2. Here we fix the octupole strength >q = 1 and the pairing strength 
equal to the energy spacing between the single-particle levels, e = G = 0.2. 

In the right column we show for comparison the spectrum for the pure pairing interaction, k<l = X3 = 0. As 
well known from the analysis of the exact solution of the pairing problem [3l| . partial seniority quantum numbers 
(or quasispins) sj are conserved, and every eigenstate is characterized, apart from the angular momentum quantum 
numbers IM, by the seniorities Sj, Sji. The states with nonzero seniority are multiple degenerate with respect to total 
spin. The states with the odd seniority Sj* have negative parity. The degenerate states with Sj — 2, Sji = show a 
clear pairing energy gap that exceeds the single-particle spacing e. 

The first four columns show spectra in the presence of multipole-multipole interactions for X3 = 1 and different 
values of The quadrupole collective state 2 + emerges within the pairing gap with energy only weakly changing 
as a function of the quadrupole constant. At relatively weak quadrupole strength, the collective low-lying 3~ state is 
also present. However, as a consequence of our limited single- particle space, the octupole and quadrupole branches 
essentially compete because, due to their opposite parity, they require very different distortion (polarization) of the 
quasiparticle vacuum. 

The situation is different in the neighboring odd system (N — 7), Fig. [21 where the same Hamiltonian parameters 
are used. In the case of the pairing only, we have two "vacuum" states with opposite parity and the same spin 17/2 
which correspond to a position of the odd particle on the j-level or on the j'-level. This asymmetry sets into action 
the mechanisms qualitatively described in Sec. 2. We see low-energy collective excitations, which, in the region of 
approximately equal strength of quadrupole and octupole branches, create parity doublets (17/2+, 17/2 - ). Their 
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Figure 2: Same as Fig. [T]but for a 7-particle system. 



minimum energy splitting, 0.15 at X2 = X3 = 1, is smaller than the single-particle spacing. 



C. Degenerate case 

In the model with j = j , the interactions are symmetric with respect to the interchange X : (J <-> j') of levels 
of opposite parity. This operation commutes with spatial inversion in an even system and anticommutes in the odd 
one, VX = {—) N TP, where N is the number of particles. For degenerate levels (e = 0), X becomes the symmetry of 
the strong Hamiltonian. But for an odd system, X-symmetry is incompatible with parity. Since the full X-symmetry 
is accidental, we still classify the nuclear eigenstates by parity II but in this case two states of opposite parity are 
exactly degenerate being related by the J-operation. Thus, in the odd system we have exact parity doublets, and all 
many-body states are pairwise degenerate. In the even system, V and X have common eigenstates with corresponding 
values ±1. 

In our definitions, both the weak interaction W and the Schiff moment operator S change parity but preserve X- 
symmetry. The transition matrix elements of S in the even nucleus occur only within the class of states with the same 
X. In the odd case we classify the stationary states (with no weak interaction) by parity so that the new symmetry 
places no constraints onto matrix elements of S. 

The spectrum of the exceptional case, e = 0, is shown in Fig. [3] for N = 6. In the pure pairing model (the right 
column), the well known spin-degenerate seniority spectrum (s = Sj + Sji) is given by 

E(s) = - , G i -r (N-s)(4j + 4-N-s), (51) 



where even values of the seniority quantum number s / increase energy through the Pauli blocking effect on 
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Figure 3: The spectra of the 6-particle system with degenerate single particle levels e = 0; G — 0.2 and >cz = 1.0. The set of 
values ^2=0.5, 1, 1.5 and 2 is used. The pure pairing spectrum, — Hz = 0, with total seniority marked is to the right. Most 
of the lowest states are labeled as J VI . The negative parity states are shown in red. 

quenched pairing. The states for non-zero multipolc-multipolc forces are labeled by quantum numbers 1TX. The 
situation with the lowest collective states inside the gap is however not significantly different from that observed for 
e ^ 0. Again as x.2 grows, we see the destructive influence of the quadrupole mode onto the octupole mode; the 
latter is strongly pushed up in energy In our analytical consideration, Sec. 2.1, the interplay of the qua drupole and 
octupole modes was constructive (and this is what is seen in global systematics and in xenon isotopes [24]]). It follows 
from eq. ^ that the mutually supportive interaction of the modes comes from the three-phonon ahnarmonicity, 
while the conventional Hamiltonian of the type (|38|) makes the coexistence of the modes less probable, at least in 
the truncated orbital space. We need to mention the crossing of the lowest 2 ++ and 2 H levels at >ti = 1.27; then 
the matrix element of the Schiff moment between the lowest states 2 H and 3~~ does not vanish, in contrast to the 
matrix element between 2 ++ and 3 . 



D. Schiff moment 

Fig. [4] shows the results for the matrix elements of the Schiff moment (in the even and in the odd systems) and of 
the weak interaction (in the odd system). In the even case (the left column, N = 6), the Schiff moment (the reduced 
matrix element between the lowest "phonon" states 2 + and 3~ shown in the upper left panel) strongly depends on 
the single-particle spacing which was partly discussed above. For all values of e, the maximum of this matrix element 
is reached in the region of approximately equal strength X2 and K3. 

The three remaining panels on the left show the excitation energies of the lowest 2 + and 3~ states for different 
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Figure 4: The system of two j — 17/2 levels is studied as a function of quadrupole interaction strength >C2 at fixed values 
X3 = 1 and G = 0.2. Three cases of single-particle splitting are considered: e = (or very small value of 0.01 to remove abrupt 
numerical discontinuities), 0.2 and 0.5. The left column is for the N = 6 particle system. Three lower left panels show the 
excitation energies of 2 + and 3~ levels for different e as labeled. The top left plot shows the reduced matrix element of the 
Schiff operator between these states. The 7-particle system is depicted in the right column plots. The splitting of the ground 
state doublet is on the right top followed below by the plots of expectation values of the Schiff operator, weak matrix element 
and their product. Three curves, solid (black), dotted (red) and dashed (green), correspond to the same set of single-particle 
spacings e — 0, 0.2 and 0.5, respectively. 



single-particle spacings e. As e and X2 increase (at fixed X3), the destructive interaction of the modes pushes the 
octupole phonon energy to high energy. However, at the reasonable value e = 0.2 that corresponds to the spectra 
shown earlier, there is a region, where X2 and ^3 are of the same order and both resulting phonon frequencies are low. 

The right column of Fig. 2] gives the results for the odd system, N — 7, as a function of w%. The upper panel 
shows the splitting of the parity doublet E{YI /2 + ) — E(17/2~). Both matrix elements, of the Schiff moment and of 
the weak interaction, show a well pronounced maximum at K2 ~ kj, ~ 1 for not very large single-particle spacing e. 
This is exactly what was the purpose of the study. Of course, the precisely degenerate case, e = 0, is unphysical (in 
order to justify perturbation theory (|25p . the actual calculation here assumed a very small spacing still exceeding the 
off-diagonal matrix element of the weak interaction). However, a noticeable effect is present in the case of e = 0.2 as 
well. 



IV. CONCLUSION 



We discussed the arguments in favor of the hypotheses of possible enhancement of the nuclear Schiff moment, and 
therefore of the effects of V , T- violating weak interactions, in nuclei with the simultaneous presence of soft collective 
modes of quadrupole and octupole nature. If this idea is correct, the pool of possible nuclei — candidates for the 
successful search of the atomic EDM — would significantly broaden. 

In our analytical model, we showed that in the odd neighbor of the even spherical nucleus with soft vibrational 
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modes the unpaired particle leads to the spontaneously broken rotational symmetry. The arising effective deformation 
can be described as a condensate of phonons. The strength of the condensate is inversely proportional to the small 
phonon frequencies in the even nucleus. This leads to the enhancement of the intrinsic Schiff moment, in analogy to 
what was known for nuclei with static quadrupole and octupole deformation. As a result, we obtain the enhancement 
of the Schiff moment in the laboratory frame that emerges due to the V ', T- violating weak interaction, whose discovery 
is the primary goal of experimental and theoretical research in this area. 

The exact numerical study of a similar two-level model confirms the existence of the region in the parameter 
space (single-particle spectrum and the strengths of pairing and multipole-multipole interactions) where the low-lying 
quadrupole and octupole modes coexist, and where matrix elements of the weak interaction and of the Schiff moment 
are indeed enhanced. The exact calculation does not introduce any intrinsic frame and fully preserves rotational 
symmetry and (in the absence of the weak interaction) parity of stationary states. Although the model studied above 
had a limited single-particle space and a small number of valence particles that precludes the appearance of very 
strong collectivity, the presence of the enhancement is visible. 

We pointed out the difference between the analytical consideration and the numerical model in the character of the 
interaction between the quadrupole and octupole modes. The numerical study, being exact in the framework of the 
two-level model, assumed the conventional two-body Hamiltonian. Then the coexistence of the modes of incompatible 
symmetry is suppressed by their destructive competition for the available supply of single-particle excitations. The 
analytical model assumed the strong anharmonic interaction between the modes that leads to their mutual support 
and the appearance of the quadrupole condensate in the state of the even nucleus with one octupole phonon. Such 
cubic anharmonicities are rather weak if considered in the next order of the random phase approximation built on the 
two-body multipole-multipole interaction. However, currently there is a broad discussion of the role of three-body 
forces, see for example 32], in the nuclear structure beyond light nuclei. The possible collective effects of three-body 
forces (in this respect it does not really matter, bare or effectively induced in nuclear matter) were briefly discussed, 
in relation to various physical phenomena in nuclear physics, in [331 ] - One of the main new effects is the emergence of 
three-phonon anharmonicity. The predicted positive correlation of quadrupole and octupole modes is well seen in the 
chain of xenon isotopes [24]. Such studies can be important in many questions of nuclear structure, independently of 
the problem connected with weak interactions. 
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